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Abstract. We perform analytical and numerical study of static spherically 
symmetric solutions in the context of Brans-Dicke-like cosmological model by Elizalde 
et al. [T] with an exponential potential. In this model the phantom regime arises 
without the appearance of any ghost degree of freedom due to the specific form of 
coupling. Eor the certain parameter ranges the model contains a regular solution which 
we interpret as a wormhole in an otherwise dS Universe. We put several bounds on the 
parameter values: w < 0, a^/|w| < 10“^, 22.7 25. The numerical solution could 

mimic the Schwarzschild one, so the original model is consistent with astrophysical 
and cosmological observational data. However differences between our solution and 
the Schwarzschild one can be quite large, so black hole candidate observations could 
probably place further limits on the (j)Q value. 


1. Introduction 

The astrophysical data (especially from the last decades) ranging from high redshift 
snrveys of snpernovae to WMAP observations indicate that onr Universe is experiencing 
an accelerating phase of expansion now [2]. A possible interpretation of this expansion 
in terms of General Relativity states that abont 70% of the total energy of onr Universe 
is attribnted to the dark energy with large and negative pressnre (cosmological constant 
is the best £t nowadays, see however [3] for a comprehensive review on dark energy 
ore [1] for a brief one). One of the simplest possibilities of a theoretical description of 
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the dark energy, beyond the cosmological constant, is to add a scalar held into the GR 
action (a scalar-tensor theory). 

Brans-Dicke (BD) scalar-tensor theory was one of the hrst attempts to modify 
gravity [5]. This theory introduces a scalar held 0 which manifests itself as gravitational 
constant variation. The action of BD theory has the form |5]: 



( 1 ) 


Through the years the theory passed multiple observational and theoretical tests. A 
variety of physical phenomena can be interpreted in terms of Brans-Dicke theory, ranging 
from the hat galaxy rotation curves [0] to inhation [7] and cosmological constant problem 
(vacuum catastrophe) [8]. It is now also well established that Brans-Dicke theory can 
be used for dark energy modelling p[ini[n]. 

The current observational bound on the BD parameter |a;| > 50000 is imposed 
in terms of the parameterized post-Newtonian (PPN) expansion of the theory |12j . 
Such a large value causes debates on the theory’s adequacy since it’s PPN expansion 
reaches general relativity (GR) in the limit u ^ oo (see [13] for detailed information 
on connection between BD and GR). Nevertheless, papers [I1]-[I7] show that most 
of the scalar-tensor theories in their cosmological evolution reach the state where the 
contribution of the scalar part to gravity virtually vanishes. So although u is large 
nowadays, it could demonstrate more interesting values (from the observational point of 
view) in the past. On the other hand, regardless of the value of u (paper [18] proposes 
|a;| > 10^^°) for |a;| < oo the theory can describe bouncing cosmology: considering the 
evolution of the Universe backwards in time, one can see that the scale factor does 
not turn to zero at any moment. Its value decreases to a minimum (bounce) and then 
starts increasing again. All the characteristic functions remain regular at the bounce 
|18] . contrary to the presence of the initial singularity in the classical GR cosmology. 
Hence the theory is of scientihc interest since large u simply provides the agreement 
with observational data and also helps avoiding the cosmological initial singularity. 

In the BD theory the possible way to take dark energy into account is to add a 
scalar held potential to the action [T^ [20] . As it was shown in m for an arbitrary 
form of the scalar held potential the Universe evolution could mimic the AGDM model 
(i.e. observational data) in BD theory. In this paper we consider spherically symmetric 
solutions in the BD-like model arising in the cosmological context in the paper by 
Elizalde et ah [1]. The action has the form: 



Here a and 0o are constants. In such a model the phantom regime arises in a BD-type 
scenario without the appearance of any ghost degree of freedom due to the specihc 
form of the coupling. It is also shown in [1] that quantum gravity ehects may prevent 
(or, at least, delay or soften) the otherwise unavoidable hnite-time future singularity 
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(Big Rip), associated with the phantom. Although this action can seem quite ad- 
hoc, it can naturally emerge in the low-energy limit of string/M-theory [211 [IIj which 
means that the cosmological phantom can emerge naturally within the multidimensional 
unihcation theories. All this makes the consideration of this model very attractive from 
the theoretical point of view. 

Elizalde et ah consider late-time cosmology with equation of state parameter value 
close to —1 (this value is favored by observations [2]). Under these circumstances 
an exact spatially-flat Friedman-Robertson-Walker (FRW) cosmology is constructed 
admitting acceleration phases for the current universe. The model agrees with the 
observational data in a wide parameters range. Such behaviour is not spoiled at the 
perturbation level since the scenario is free of perturbative instabilities |T] . All the above 
arguments make this model perspective for further versatile research. 

Negative values of uj are often ruled out because of the wrong sign before the kinetic 
term in the action leading to a ghost. Nevertheless a plenty of researchers consider a; < 0 
for the following reasons. The latest observational data indicate that phantom nature of 
dark energy is more likely |2]. Combining WMAP-|-eCMB-|-BAO-|-HO-|-SNe data yields 
to wdeo = —1.17(-|-0.13—0.12) for the flat Universe at the signihcance level of 68 %. Such 
a trend also occurs for a non-flat Universe model and for different dataset combinations. 
Finally, consideration of a scalar held as an effective description for the theory with 
positive dehned energy could eliminate the quantum contradictions [ 22 l [23]. All these 
arguments make the consideration of BD theory with ca < 0 reasonable, since it leads 
to a phantom cosmology. Modern observational bounds limit the absolute value of uj. 
Modelling the gravitational collapse does not rule out negative u values as well [2l] . So 
the range a; < 0 is of great interest from the phenomenological point of view, since it 
can provide no-ghost but phantom cosmology in agreement with modern observations. 

In this paper we explore static spherically symmetric solutions of (|2]). Cosmology 
studies scalar held time dependence 0 = 0 (t), this approximation stands for cosmological 
homogeneity. Understanding the local ehects 0 = 0(r) yields to a picture of 
held behaviour at diherent scales and checks for model consistency. Exploring the 
cosmological model from another point of view is a step towards a better understanding 
of dark energy. 

The most well-known static spherically symmetric solution in BD was obtained by 
Campanelli and Lousto [25]. This solution was independently rediscovered by Agnese 
and La Camera m and correctly reinterpreted as a naked singularity for 7 < 1 and 
a wormhole for 7 > 1 (here 7 is a PPN parameter). Scalar held plays the role of the 
exotic matter at the wormhole throat and ensures it’s traversability not only for a; < 0 
but also for large positive uj [ 2 ^. Astrophysical properties were studied by Alexeyev et 
al. in [28] and found to be in agreement with modern observations. 

The purpose of this paper is to explore the properties of the static spherically 
symmetric solutions in the framework of Elizalde et al. To do this we analytically 
dehne the parametrization for the pure Agnese and La Camera solution and obtain a 
new analytical solution for the specihc metric ansatz, which can be called a “stealth 
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Schwarzschild solution” (analogously to [21] )• The presence of the scalar held potential 
modihes the held equations, so new types of solutions can emerge. We perform a 
numerical exploration of this solutions and discuss their properties. Elizalde et ah 
construct a model in such a way that the ghosts can be absent, so it is interesting to 
check whether regular local solutions can appear in the ghost-free sector of the model. 

Spherically symmetric solutions for similar actions are present in the literature. An 
asymptotically Lifshitz black hole solution for a power-law potential in BD setting was 
found in 1301 but the allowed u parameter range appears to be too narrow to ht the 
observations. Another wormhole solution for the potential of the form V{(j)) = A0 is 
given in [31]. The latter one is obtained for cj) = this wormhole is non-static, its 
throat radius will increase as the cosmic time increases. For the best of our knowledge 
these are the only solutions in the BD setting with a scalar hied potential. The paper 
[32] considers quantum wormhole solutions for the BD gravity with A in the absence 
of matter perturbatively, from the held theory point of view. Although the framework 
of [32] is similar to our and can be obtained via V{(j)) = 2A, the results of [32] cannot 
be applied to astrophysics. Wormhole in solutions the GR setting with exotic matter 
and A representing similar behaviour were found in [3^ IM] . The cosmological constant 
ahects the wormhole and its properties, the wormhole geometry will be dominated by 
the de Sitter space far from the throat. Our solution is new and dihers from the ones 
of [33l IM] since the role of the exotic matter is played by the scalar held. 

The outline of our paper is the following. In Section |2] we study the generic 
properties of discussed BD-like solutions. Section [3] is devoted to specihc analytical 
solutions, in Section 0] we obtain some new limitations on model parameters. Section 
|5] is devoted to ISL modihcations, in Section |6] we show the results of numerical study. 
Section [6] contains discussion and conclusions. 


2. General properties 


The corresponding Einstein and Klein-Gordon equations for the action ([2|) have the 
following form: 

an + uad^(f)d^(t) - + 2cnn0 = 0 , (3) 

G^,=ul - lg,,da(j)dy] + e-“^ - g^^Ue^^) - lg,,V{(j)). (4) 


Hawking’s theorem on BD theory [35] is usually taken to mean that BD static 
spherically symmetric solutions are exactly the same as those of GR. We want to note 
that this is an overstatement. For example, the black hole solution of Gampanelli and 
Lousto [25] is spherically symmetric and static but does not match the Schwarzschild 
one. The proof of Hawking’s theorem relies on the equivalence between Jordan and 
Einstein frames and the fact that the corresponding conformal transformation is well- 
dehned. As it is also pointed out in [3^ the advantage of using the Einstein frame is 
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that in the Einstein frame BD scalar has canonical kinetic energy density and obeys the 
weak (WEC) and null (NEC) energy conditions. If the conformal transformation to the 
Einstein frame becomes ill-dehned at the horizon (as it occurs in [25]) then its variables 
dEfiu, <pE (subscript E stands for Einstein frame) cannot be used on such surface [36|- The 
last statement invalidates the proof of Hawking’s theorem. One could use, of course, the 
Jordan frame instead of the Einstein one. In this case the scalar 0 violates the WEC and 
NEC because its stress-energy tensor has a non-canonical structure containing second 
derivatives of 0 instead of being quadratic in the hrst derivatives. 

So, if our Einstein frame scalar does not obey the WEC or the Jordan frame scalar 
diverges at the horizon then the Hawking theorem does not apply and the solution is 
not forced to be Schwarzschild one. 

One can perform the transformation to the Einstein frame g^uE^"^ = = dE^iu- 


Sr,= 


IBtt 


\/-gE 


ds% = e'^^dt% — — r^^{d9^ -|- sin^ 6d(f^) , 


TZe 

2 


ujd^(pd^(t) - V ( 0 ) 


CJ 


= (3a V2 + Cj), H(0) = Ve-^^ = 


(5) 

( 6 ) 
(7) 


It follows from the expression for cj ([7]) that even if cj is negative, in the case when 


“io? J- ClJ > 0 


( 8 ) 


the effective kinetic energy of 0 becomes positive, similarly to the usual scalar held, and 
therefore the ghost does not appear, although the cosmology is phantom [1]. Introducing 
the four-velocity 0, 0, O) we obtain the WEC in the form 

a;0'^/0^e“'^'® — V{(j))> 0, (9) 


which is generally satished everywhere if Veo < 0 and 3a^ -|- 2a; > 0. 

Via redehnition <h = e“'^ one can obtain the theory in the form of the standard BD 
one plus potential: 







a^ <h 






( 10 ) 


Having in mind that we are now working with the “usual” BD we can directly apply 
it’s well known weak held limit results to our model. We obtain that G ~ 
where G is the ehective gravitational constant. 


3. Analytical solutions 

Let us rewrite the held equations as in [25]. The Klein-Gordon equation ([3]) can be 
rewritten as 

aR - ( « + — 1 + —= 0 . 

\ 00 / a 


( 11 ) 
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Multiplying (jl]) by g^’’' we obtain 

R - - 2Voe^/^° + = 0 . 

Combining equations (IT^ and (ITTll we get 

^ + Sa") - Voe^^^° (a-^]e^‘^ = 0 . 


( 12 ) 


(13) 


Substituting (IT^ and flT^ into the field equation (jl]) we arrive at the following equation 




D = Vn 


a 


a 


2u + 3a2 


(14) 

(15) 


Further we obtain specific analytic solutions for fixed parameter relations or specific 
metric structure. Of course, the general analytical solution remains unknown. In section 
[ 6 ] we investigate the general solution numerically. 


3.1. The case D = 0 
When H = 0, or, equivalently 


00 — 


a 


(16) 


2(0)2 + (u)’ 

one has the well-known solution of Agnese and La Camera (ALC) [26], since the field 
equations match in this case. Within the framework of Elizalse et ah this solution looks 
like 


ds^ = A{r)^df 


n — m 


-r o -In 

^ 2 a 2 

V = M\/(l + 7)/2, 



A(r)i-V(dr+ sin^ 0 V), 


2rj 

A{r) = 1 - 

r 


m 


ui 

~~2 


m 


2n 


a 


n — m 


G = e- 


= y2/(i + 7), 

, 2 


n 


1 + 7 ’ 


( 17 ) 

(18) 

= 772/(1 + 7 ), (19) 

( 20 ) 


with 7 being the PPN parameter, 7 plays a role of the mass parameter, G is the effective 
gravitational constant and F can be fixed from the Newtonian limit. The parameter a 
comes from the coupling of Elizalde et ah action ([2]). This solution corresponds to a 
naked singularity if 7 < 1, to a Schwarzschild black hole for 7 = 1 and to a wormhole 
if 7 > 1. The scalar field diverges at the horizon and this is the reason of the Hawking 
theorem violation. The factor 1^(1 -|- 7)/2 is usually absorbed in the mass definition, 
which results in the correction of order of 10 “® for I 7 — 1 | ~ 10 “^. 

It is worth noticing that the condition D = 0 preserves the potential term in the 
action but effectively removes it from the Einstein equations. If we suppose that the 
scalar field potential acts as the cosmological constant, we can arrive at a conclusion 
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that the accelerated expansion of the Universe does not contribute to local dynamics. 
This effect is indifferent to the value of the cosmological constant. 

To estimate the constant F we must set up the Newtonian limit, in other words, 
get the model linearised. This work is shown in the Appendix. Comparing (l55|) with 
flTSl) for large distances we arrive at 

2a 

1 2u + 4q!^\ n — m 

For the power-law generalization of Schwarzschild flTTl) the equations flTdl) can only be 
solved for m = n which leads to D = 0. So there are no other power-law solutions for 
the ansatz flTTl) except ALC here. 



1 2a; -|- 4a^ 
^2u + 3a^ 


F = 


3.2. Schwarzschild-like solution 

The other solution can be obtained if we suppose that we measure exactly Schwarzschild 
metric 

_ gA ( r )^^2 _ ^ ^-^^2 ^^^ 2 ^ ^ ^ 22 ) 

z/= ln(l - r/j/r) A = - ln(l - r/i/r). (23) 

This leads to the solution 

a 

4> = \r + Ci\ + C 2 ,r ^ -Ci (24) 

a^ + a;2 

Having in mind the redehnition flTOl) we see that G is a decreasing function of r 

G ^ {r + Cl) 1+“^/“^. (25) 


This solution can be called a “stealth Schwarzschild solution”, analogously to the 
solution of [29] (obtained within the Horndeski/Galileon framework) representing a 
Schwarzshild metric with non-trivial scalar held. 

To estimate the constants on fl25l) we must use the Newtonian limit, in other words, 
get the model linearised. The scalar (j) increases indehnitely upon r, so the theory can 
only be formally linearised as 0 = 0c + ^ in fhe region outside the object ^ < 1 where 
the Taylor expansion of logarithm exists. The constant Ci represents then the length 
scale up to which the solution admits Newtonian behavior. 

The obvious condition preserving (to a certain extent) the observational adequacy 
of the solution is: |a;| 3> |q;|. It makes the power of r small and ensures that ip does not 
change abruptly space and varies less and less with the distance. This condition can 
also be agreed with fl27|l . The space-time is asymptotically flat, but the gravitational 
constant vanishes at the inhnity. 

The discovery of the discussed solution under the assumption that it realizes in 
nature, seems tricky, since the metric is exactly the Schwarzschild one. It only manifests 
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itself as a gravitational constant lowering itself with distance. One could naively suggest 
that for astrophysical scales this effect would simply hide itself into mass dehnition, since 
accurate mass measurements are done mostly for several binary systems due to their 
orbital motion [37]. Therefore a lower but very slowly changing gravitational constant 
value would just be absorbed into calculated masses: the less would be the distance, 
the larger the object masses would seem to us. This situation contradicts to the dark 
matter effect: if we suppose it to be dynamical and not caused by some real-life matter. 


4. Parameter estimations 


From the paper by Elizalde et ah [T] we can get the cosmologically imposed condition, 
providing the accelerated expansion of the Universe: 



4 

3a2 + 2n;' 


(26) 


For |ci;| 3> \a‘^ \ and u: is negative this condition is always satished and is also in agreement 
with (ITbl) . This means that ALC solution can represent a wormhole embedded in an 
expanding Universe, but not affected by this expansion at all. 

The other possible parameter estimation comes from the PPN expansion. The 
observational bound is [12] 7 — 1 = (2.1 ± 2.3) x 10“®. Following [3H] we derive the 
PPN-parameters 7 and /9: 


7-1 


2a^ -f a;’ 


/ 3 - 1 = 0 


When I a; I S> |a^|, taking into account I 7 — 1| < 10 ® we obtain 


(27) 


a^/|a;| < 10 


(28) 


This requirement disagrees with no-ghost condition ([H]) for negative uj. In fact the 
inequality above leaves no room for the no-ghost sector of the model. This result is a 
bit disappointing, but is not critical in view of the reasons regarding negative u: given 
in the introduction. 


5. Inverse-square law modifications 


The power-law modihcation of the inverse-square law (ISL) in the general form reads 


5J(r) 


GoM 


(W" 

1 I “ j 


r 


(29) 
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Expanding the ALC (fT7|) - fl2n]) solntion gQQ in binomial series and going to physical nnits 
we get 


QJ 


T]C 


1 + 7 r 


1 + 



GqM 

r 


7 ] 

1 + 02 — ) 
r 


(30) 


02 — 



V = 


'I + 7MG0 


(31) 


for = 45 . We are interested in the case when 7 > 1, so 02 > 0. It is worth noticing 
that the correction is mass-dependent. Papers [iQl [HI |12] claim the following Icr limits 
for power-law corrections 


N = 2 ^ < 1.3 x 10 ®[m]. 


(32) 


For the conditions of the experiment |1T] (M ^ 10kg) we estimate 


027 


MGo 




e X 10"^®. 


(33) 


Together with fl32l) this gives ns e < 10^®. For the PPN parameter I 7 — 1| < 2.5 x 10“® 
we can estimate e < 2.5 x 10“® which is a snfhciently stronger bonnd. This means that 
one needs to increase the experiment accnracy by 16 orders to approach the PPN-based 
limit. The satellite motion in the Earth gravitational held or the planet motion aronnd 
the Snn conld represent a more interesting case whereas Mq iW^kg, Me ~ iW^kg. 

Following [13] we estimate additional relative freqnency shift (5///)^ for the 
potential above and make bonnds nsing the freqnency measnrement accnracy for Galileo 
Navigation Satellite System (GNSS)@. 

^ ( I i t KGl f 1 i t 

V/JA (Rti + hf) c* (R^ + hf) 

Gonsidering the accnracy of the freqnency measnrement of the Galileo constellation 
= 10 “^^ analogonsly to fl3^ we obtain 

<5///^5exl0-'3<£j,. (35) 

f The interesting paper [39] is worth mentioning here. The authors calculate advances in the perihelia of 
planets in the solar system by treating the power-law correction as a small disturbance and then connect 
them with the data of INPOPlOa (IMCCE, France) and EPM2011 (lAA RAS, Russia) ephemerides(by 
means of minimizing the sum of squared deviations). They take the uncertainty in the Sun’s quadrupole 
moment into account and estimate it along with the parameters of the power-law correction. The result 
is N = 0.605 for the exponent. However, from EPM2011, they find that, although it yields N = 3.001, 
the estimated uncertainty in the Sun’s quadrupole moment is much larger than the value of ±10% 
given by current observations. So, no any certain conclusions on the most plausible value of N can be 
done. 

§ h = 23.222 X lO^/cTO, Me = 5, 972 x lO^^kg, R® = 6371A:to, Go = 6.67384 x lO-^^m^kg-N-2 
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The resulting bound appears to be as e < 2. It is also much weaker than the PPN-based 
one. This result seems strange, since the experimental technique which provides the 
PPN and the GNSS bounds is similar. This peculiarity requires additional research. 


6. Numerical results 


Now we intend to analyse the general solution. For numerical calculations we use the 
following ansatz 


— -^dr^ — R{rYdVl? 


(36) 


Without loss of generality we set a = 1 since this is equivalent to the radial coordinate 
transformation. From the equation flTT|) we obtain the system of differential equations, 
which can be resolved with respect to highest derivatives as 




A'Tr) = 


B"(r) = 


DR{r)e — 2A'(r)i?'(r) 

mr) 




-2R{r)A'{r)R'{r) — 2A(r)R'{r)^ + DR{r)^e + 2 


2A(r)i^(r) 


(37) 

(38) 




„ + O') 0'(r)^ 2 (A(r)i?'(r)^ — 1) DR{r)e‘I’o — 2A'(r)i?'(r) 

^ ^ ^ aA(r)i?(r)2 aA(r)i?(r) ' 

The most interesting scales from the phenomenological point of view are the 
astrophysical ones, such as neutron stars binary systems. Hence we will consider a 
mass of order of the solar one. 

Under the requirement for the potential term in the action to transform to a 
cosmological constant in the limit r ^ cx), 0 ^ 0oo we obtain that 


Uo = 2Ae“'^“/'^°, 


0OO 


1 

— In 

a 


1 20; + 4q;^ 
^2a; + 3a2 ’ 


(40) 


where A is the cosmological constant. Asymptotic value of the scalar held is derived 
from (155|) . 0o is a free parameter. The initial values for the numerical integration are 
taken from the ALC solution. Summarizing the known bounds we can state that o; < 0 
(no naked singularity in ALC), < 10“^|o;| (PPN bound). For the sake of illustration 
we take u = —10000, a = 0.1, M = mJJJ]. 

Since the values of oj and a under consideration provide WEC violation, and 
requiring the correspondence with Agnese and La Camera solution for small potential 
inhuence we claim that the solution represents a wormhole embedded in an otherwise 
de Sitter universe (see hg. [1]). Near the throat, the wormhole geometry is dominant 
while as the distance increases the structure of dS reveals itself. 

The space-time approaches the Schwarzschild metric when one continues increasing 
00 (this corresponds to reducing the potential influence) and differs slightly from 

II c = 2.998 X 10® m/sec, Mq = 1.988435 x 10®® g 
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A 



Figure 1. The solution represents a wormhole embedded in an otherwise de Sitter 
universe ((/jq = 24.5). Near the throat, the wormhole geometry is dominant while as 
the distance increases the structure of dS reveals itself. 


Schwarzshild for 0o > 25. For 25 < 0o ^ 22.7 the solution represents a dS- 
wormhole. Lowering 0o further (strengthen the potential influence)we see that space- 
time approaches dS. For 0o ^ 22.7 the wormhole throat does not form, suppressed 
by the potential term influence (dS case, see £g. |2]). De Sitter solution describes the 
expanding universe. Wormhole-like solution describes an astrophysical object. Thus 
the solution can describe both cosmological and astrophysical scales. 

A 

8 - \ 

\ 

\ 

7 - 
6 - 
5 - 
4 - 
3 - 
2 - 

0 


\ , 

\ - 
\ ■ 
\ 


100000 200000 



Figure 2. Figure represents the enumeration of do values. Approximately at do ~ 22.7 
the solution turns from de Sitter to a wormhole-like. De Sitter solution describes the 
expanding universe. Wormhole-like solution describes an astrophysical object. Thus 
the solution can describe both cosmological and astrophysical scales. 

Figures [3] and m display that the given numerical solution is a wormhole analogue 
of the Schwarzschild-dS solution. Near the throat, the wormhole geometry is dominant 
(fig. [3]) while as the distance increases the structure of dS reveals itself (fig. Ilj). The 
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threshold on 0o between dS and wormhole-like geometry appears near 0o = 22.7 (See 
Fig. E]). 



Figure 3. Figure represents the enumeration of do values. The threshold on do 
between dS and wormhole-like geometry appears near do = 22.7. 


A 

12 - 

10 - 

8 - 

6 

4- 

_ ' ■ ■ • ■ 27 

_ ■ — ^ 

-2.x 10^ -1.x 10^ 0 


% 

- 23.2 

- 23.525 

- 24 

24.5 



6 6 I* 

l.xlO 2.x 10 


Figure 4. The solution represents a wormhole embedded in an otherwise de Sitter 
universe (do = 24.5). Near the throat, the wormhole geometry is dominant while as 
the distance increases the structure of dS reveals itself. 

Solutions representing similar behaviour were found in [331 l3l] . The papers [331IM] 
consider the GR setting with exotic matter and A. The cosmological constant affects 
the wormhole and its properties, the geometry will be dominated by the wormhole near 
the throat and by the the de Sitter space far from the throat, as it is for our solution. 
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Our result are in agreement with [Ml El], however the role of exotic matter is played in 
our case by the scalar held. 

Resetting the metric in a form 

dr‘^ 

ds^ = - . - r'^dVt^ , (41) 

1 — b[r)/r 

we can plot the shape function h{r) and the embedded surface z{r) |M] 

dz ( r V'''" 

To be a solution of a wormhole, the geometry has to have a minimum radius, r = 
6(r) = ro, i.e the throat, at which the embedded surface is vertical dz/dr —>■ oo 
(see hgure [5]). For a wormhole to be traversable one must demand the absence of 


z(r) 



Figure 5. The embedded surface is vertical dz/dr — >■ oo near the wormhole throat for 
bo = 24.5. 

horizons identihed as surfaces with —)■ 0, so that the function P{r) has to be hnite 

everywhere. By the proper choice o R'{r) it is possible to avoid such a divergence at the 
throat. The numerical calculation shows that for 0o = 24.5 the throat radius is around 
r = 5927m, this is dehnitely large enough for a microscopical object to pass through. 
The corresponding ALC throat is at r = 3014m, so the potential influence is signihcant. 

The difference between the Schwarzschild solution and numerical solution at the 
distance around the innermost stable orbit r = 6 y is 

Nnum/^schu, ~ 1.031 0o = 24.5, (43) 

^num/^schw ~ 1.471 (pQ = 23.5. (44) 

These values appear to be quite large, so black hole observations could probably place 
some limits on 00 - There are observational indications that such a deviation might in 
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fact correspond to reality: the paper [H] claims that the apparent size of the source at 
the galactic centre is less than the expected apparent size of the event horizon of the 
presumed Sgr A black hole. If this is not due to modelling issues, this could be thought 
as the indication of GR violation. To explore the properties of general solution further 
in detail one requires the analytical form of the solution and we treat this as a task for 
the future. 

We need to distinguish between the black hole, wormhole and naked singularity 
cases. We check different parameter combinations and, compare Kretschmann curvature 
invariant (see Fig. |6]) and metric function A (1361) values. The analysis of multiple 
plots shows that only the case a; < 0, 22.7 ^ 0o ^ 25 represents a regular solution: 
Kretschmann scalar is hnite at r > 0. Therefore this case can represent a traversable 
wormhole. 

All other parameter ranges seem to represent a naked singularity for different values 
of r (see Fig. [6] also). The cosmic censorship conjecture states that naked singularities 
do not emerge in nature. In terms of wormhole research naked singularities can be 
interpreted as not traversable ones, so they are of a limited interest. However nowadays 
the cosmic censorship conjecture remains unproven and many researchers suggest naked 
singularities as astrophysical object candidates [l5l HG]. In this sense it would be 
interesting to check for the astrophysical properties of the solution in corresponding 
parameter regions. We live it for the next step of our investigation. 

K 

-15 

3. X 10 - 


2. X 10 - 

-15 

1. X 10 - 



5000 10000 

Figure 6. The Kretschmann scalar curvature invariant and the metric function A for 
different parameter values. Only the case w < 0, 22.7 ^ bo ^ 25 represents a regular 
solution. All other parameter ranges represent a naked singularity for different values 
of r. 
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7. Discussion and conclusions 

We explored analytical and nnmerical properties of static spherically symmetric 
solntions in the context of BD-like cosmological model. For a certain parameter range 
O' < 0, 22.7 ^ 00 ^ 25 the model contains a regnlar wormhole solntion. The space- 
time geometry is analogonsly to Schwarzschild-dS one, representing a wormhole in the 
dS Universe. The inflnence of the potential term on the wormhole is shown to be be 
signihcant. The properties of this newly obtained nnmerical wormhole solntion reqnire 
fnrther versatile research which is a task for onr next paper. 

As it was mentioned in 071 the Jordan and Einstein frames for the theory are not 
eqnivalent from the cosmological point of view. It was also argued in [18] (for BD thery 
without the scalar held potential) that wormhole solution in the Einstein frame are 
absent unless Co = 3a‘^/2 -f-ca < 0 (in our notation). As it follows from the equation ([9]) 
the WEC can be violated due to the potential alone and one does not need to violate 
it by setting a) < 0 by “brute force” [IH]. So, unlike the case of |1H] wormholes are in 
principle possible in both frames, although the detailed discussion is only possible with 
the analytical solution in hand. However the observational bond fl28ll implies that the 
WEC is violated. 

From the beginning the discussed model had four free parameters. Original paper 
P provides one useful limitation on 0o (ESj). We estimate the parameter Vq using 
cosmological constant and currently measured gravitational one values fITOD reducing 
the number of free parameters by one. We also constraint a and u parameters by fl28ll . 
When considered together equations a; < 0, 22.7 0o ^ 25, fl26|) . fl28|) and fHOjl put 
signihcant constraints on the parameter values in compared with p. 

We demonstrate that when the equations (|26il . (|28ll and (l4n|l are taken into 
account the model gives rise to a potentially traversable wormhole or a naked 
singularity. Moreover, since the numerical solutions can mimic the Schwarzschild one, 
the original model P is consistent with the astrophysical observational data. The last 
statement proofs that the model could be successfully applied both on cosmological 
and astrophysical scales. However diherences between the numerical solution and the 
Schwarzschild one can be quite large (an example is the factor 1.471 in fl4^ - (144)) ). so 
black hole candidate observations could probably place additional limits on 0o in a 
case if the analytical solution would be obtained. The question is whether the proposed 
parameters can provide deviations visible for the Event Horizon Telescope (for example) 
and it requires further investigations. The numerical plotting does not allow one to 
explore all possible parameter combinations, so there ape probably other parameter 
rangers, where regular solutions exist, and our result represent just an example of such 
solutions. 

This PPN-based requirement fl28ll disagrees with no-ghost condition (0) for negative 
Lo and in fact leaves no room for the no-ghost sector of the model. This result is a bit 
disappointing, but is not critical because of the reasons regarding negative u given in 
the Introduction. 


Wormholes and Naked Singularities in Brans-Dicke eosmology 


16 


We also derive the power-law correction to the gravitational potential for the 
solution of Agnese & La Camera [26] and check whether it can be observed in the ground- 
based experiments ore via measuring relative frequency shift for Galileo Navigation 
Satellite System (GNSS). We conclude that one needs to increase the experiment 
accuracy to approach the largest possible correction value. 

We hnd out that ALC solution can represent the wormhole embedded in an 
expanding Universe, but not affected by this expansion. 

In summary, we conclude that BD theory being an effective limit of f(R) gravity 
|49] . describes a wide range of observed phenomena (ranging from cosmological to 
astrophysical scales), hts the observational data quite well, removes cosmological initial 
singularity and provides a scale factor bounce in the early Universe. So taking into 
account all the other results on BD study we would like to state that it is one of the 
hnest candidates for the foundation stone of modern cosmology and extended gravity. 
Of course the scalar-tensor theories are not free of problems, especially when they are 
directly considered as dark energy candidates. Nevertheless the attention to phantom 
models is driven by the lack of a good theoretical understanding of the present universe 
in terms of more usual theories. The phenomenology which emerges from these models is 
rich and pertinently at times, so this models deserve to be investigated comprehensively. 
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Appendix I 

Equation fl45|) in the presence of matter transforms into 

1 f2uj + 2a^\ 
T^u - -g,. y2u + 3a^ ) 

In the weak held limit we can write down 


R^i/ ujd^(l)dp(l) “t“ 


1 ,,, Svr 

—— -f -F —- 


(45) 


0 = 0c + ^(r), 




< 1 , 






□e“'^ ^ g^^ = N (46) 
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where 7]^^ is Minkowski metric and are perturbation. Up to the first expansion order 
one obtains: 


+ VuuVi + 87re-“'^' 




Ui = 


C = 


'2u} + 2W-'' 

2bj + 3a2 


On the other hand 


9p(l) _ riu \ h'^ — 

~ “T 


(47) 

(48) 

(49) 


To simplify flT7|) we impose the following conditions, using the weak held version of 
harmonic coordinate conditions in a case of the constant scalar held: 


-h^ _ 

hL,^,U 


-2a'\/ 


The equation fITTI) rewrites as 


1 


+ STre-"-^^ 


T, 


fll' 




(50) 


(51) 


Let us take the space-time and the scalar held generated by a static point-mass M. 
Therefore the corresponding stress-energy tensor is given by 


T; = diag{M5{r), 0, 0, 0), = M5{r) 


Therefore the metric perturbation reads as: 

, uj 2o? 

□hoo = 2Ui + ^ . Mb{r) 


2u + 3a2 


2G'oM Uir2 

hoo =- 1 - 5 —, 

r 6 

^ 2a; + 4a2 

Gn = 


2a; + 3a2- 


Finally we obtain the weak held limit for the considering metric in the form 


(52) 

(53) 

(54) 

(55) 

(56) 


2GoM Uir^ 

5'oo = 1-^-^ 

r 3 


(57) 


This result agrees with the fact that one can consider the model [T] in the form of 
the standard BD (ITOil where Go plays the role of the ehective gravitational constant. 
Such solution reduces to the usual BD one when a = 1, (^ = e''^, Vq = 0. After compar¬ 
ing it with the de Sitter-Schwarzschild metric we conclude that Vi plays a role of the 
cosmological constant. 
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